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Electron-impact total ionization cross sec- 
tions for the CH radical and C 2 H 2 
(acetylene) have been calculated using the 
Binary-Encounter-Bethe (BEB) model. 
The BEB model combines the Mott cross 
section and the asymptotic form of the 
Bethe theory, and has been shown to gener- 
ate reliable ionization cross sections for a 
large variety of molecules. The BEB cross 
sections for CH and C 2 H 2 are in good 
agreement with the available experimental 
data from ionization thresholds to hun- 
dreds of eV in incident energies. 
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1. Introduction 



The Binary-Encounter-Bethe (BEB) model [1] has 
been successfully used to calculate the total ionization 
cross sections of a wide range of molecules [2-4] in- 
cluding H 2 , 3 , and SF 6 . The BEB model combines a 
modified form of the Mott cross section with the asymp- 
totic form of the Bethe theory (i.e., high incident en- 
ergy) for electron-impact ionization of an atom. 

These two theories are combined in such a way that 
not only the ionization cross section has the proper 
asymptotic form but so does the stopping cross section, 
which is the product of the cross section and the energy 



transfer from the incident particle to the target [1]. Al- 
though the theoretical model in its differential form 
requires the knowledge of continuum dipole oscillator 
strengths df/dE, we assumed a simple form for df/dE to 
facilitate the integration over the ejected electron en- 
ergy. The resulting expression for the ionization cross 
section per molecular orbital, which we refer to as the 
BEB model, has a remarkably simple form [1]: 

S [lnf/. 1 \ , 1 lnf 1 ,„ 
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where t=T/B, u=U/B, S = 4Tta 2 NR 2 IB 2 , a is the 
Bohr radius (= 0.5292 A), R is the rydberg energy 
(= 1 3.6057 eV), Tis the incident electron energy, and N, 
B , and U are the electron occupation number, the bind- 
ing energy, and the average kinetic energy of the orbital, 
respectively. 

In this article, we compare the BEB cross sections for 
the CH radical in its ground state 2 I1 and C 2 H 2 also in the 
ground state ! 2g + with available experimental data. Both 
molecules are used in modeling the behavior of edge 
plasmas in the divertor region of a tokamak [5]. 



2. Results and Discussion 

The molecular constants B and U were calculated 
using the 6-31 1+G(d,p) basis set in the GAMESS code 
[6]. Although one can in principle use theoretical values 
for all ZTs, theoretical values for the lowest ionization 
potential (IP) are usually inaccurate. Therefore, we used 
experimental values [7] only for the lowest IP to match 
the observed thresholds. Theoretical results and com- 
parisons to experimental data presented below are also 
available on a NIST web page for electron-impact ion- 
ization cross sections at: 

http://physics.nist.gov/PhysRefData/Ionization/Xsection.html 

Note that the above address is case sensitive. 



2.1 CH ( 2 n) 

The molecular orbital constants B, U, and N for the 
CH radical are listed in Table 1 , and the BEB model is 
compared in Fig. 1 to the only available experimental 
data by Tarnovsky et al. [8]. Although the experiment 
was performed using deuterium instead of hydrogen, 
i.e., with CD, our theory is insensitive to the use of the 
isotope. As Tarnovsky et al. [8] have pointed out, their 
experiment did not detect fast deuterons produced, 
which may partly account for the difference between 
theory and experiment near the peak. Otherwise, theory 
and experiment are in excellent agreement, particularly 
for T < 30 eV. This is the region of T for which most 
other theories have difficulty in reproducing observed 
cross sections. The BEB model produces reliable cross 
sections for many molecules at low T (~ 50 eV or less) 

Table 1. Molecular orbital (MO) constants for CH ( 2 II): B = binding 
energy, U = kinetic energy, N = occupation number. The B value 
marked by * is an experimental value 
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B(eV) 


£/(eV) 
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Fig. 1. Comparison of the BEB cross section to experiment for CH. Solid curve, BEB cross 
section; circles, experimental data by Tarnovsky et al. (for CD) [8]. The ordinate is the 
ionization cross section, and the abscissa is the incident energy. 
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2.2 Acetylene, C 2 H 2 (%+) 

The molecular orbital constants for acetylene are 
listed in Table 2, and the BEB model is compared in Fig. 
2 to the available experimental data [9-1 1]. The experi- 
mental data by Zheng and Srivastava [9] and those by 
Tate and Smith [10] agree well with each other but they 
are approximately 15 % higher than the theory at the 
peak. On the other hand, the data by Gaudin and Hage- 
mann [11] are in excellent agreement with the theory 
from T= 100 eV to 2 keV. Since the BEB model pro- 
vides a convenient analytic form for the ionization cross 
section for the entire range of T shown in Fig. 2, and the 
data by Tate and Smith agree well with the theory for 
T < 30 eV, we are confident that the BEB cross sections 
will be appropriate and useful for plasma modeling ap- 
plications. 

Table 2. Molecular orbital (MO) constants for C 2 H 2 
( ! Xg + ): B = binding energy, U = kinetic energy, N = occupation num- 
ber. The B value marked by * is an experimental value 



MO 


B(eV) 


£/(eV) 


AT 


lCTg 


305.62 


435.15 


2 


lo- u 


305.50 


436.31 
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2o- g 


28.18 


49.60 
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2ct u 


20.80 


32.79 


2 


3o- g 


18.55 


33.64 
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1tt u 


11.4* 


28.99 
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Fig. 2. Comparison of BEB cross section to experiments for C 2 H 2 . Solid curve, BEB cross 
section; circles, experimental data by Zheng and Srivastava [9]; triangles, data by Tate and 
Smith [10]; diamonds, data by Gaudin and Hagemann [11]. The ordinate is the ionization 
cross section, and the abscissa is the incident energy. 
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